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FOREWORD 


This  report  Is  Part  II  of  the  Final  Technical  Report  of  a study 
being  conducted  by  the  Electrical  Engineering  Department  under  the 
auspices  of  the  Engineering  Experiment  Station  of  Auburn  University. 
This  technical  report  Is  submitted  toward  fulfillment  of  the  require- 
ments prescribed  In  AFOSR  Contract  F44620-76-C-0054. 
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I.  INTRODUCTION 


Numerous  Investigators  have  studied  the  o1ectvot.\njnet1c  Interactions 
of  thin  cylinders  over  a perfectly  conducting  ground  plane  Illuminated 
by  a distant  source  [1-6].  There  have  also  been  numerous  studies  of 
currents  Induced  on  Infinite  cylinders  or  thin  wires  near  a finitely 
conducting  ground  plane  [7-9].  In  this  study  a perfectly  conducting 
cylinder  of  finite  length  near  a finitely  conducting  ground  plane  Is 
treated  by  the  Singularity  Expansion  Method  (SEM). 

A Pocklington  type  Integro-differentlal  equation  Is  formulated  for 
the  current  Induced  on  the  thin  cylinder  and  Its  Image  where  the  current 
on  the  Image  Is  related  to  the  current  on  the  object  by  the  ratio  of  the 
complex  reflection  coefficient  for  the  appropriate  angle  of  Incidence  and 
polarization  Involved.  In  a recent  study  by  Sarkar  and  Strait  [10]  it 
was  shown  that  the  above  method,  termed  “reflection  method"  gave  results 
In  the  real  frequency  domain  within  lOt  of  the  exact  Somerfeld  formulation 
for  a horizontal  electric  dipole  as  long  as  the  dipole  was  at  least 
(0.25/ ^/~r}x  from  the  ground  plane.  The  advantages  of  the  reflection 
method  are  speed  of  computation  and  Interpretation  of  results. 

Using  this  “reflection  method"  In  the  Integro-differentlal  equation, 
this  study  determines  some  of  the  effects  of  ground  conductivity  and 
permittivity  upon  the  natural  resonances  and  natural  modes  of  a thin 
finite  length  cylindrical  scatteror  above  a lossy  half  sp:ce.  These 
effects  are  presented  In  terms  of  a family  of  trajectories  In  the  complex 
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frequency  plene  as  a function  of  the  electrical  and  geometrical  parameters 
^ which  define  the  problem. 

1 

1 

I 

f 

i 

X. 

I 

I 

I 

I 

I 

I 

i 

li 

D 

B 

il 


i 


a 

t 


I.  THEORY 


DERIVATION  OF  INTEGRAL  EQUATIONS 

In  Figure  2-1  an  Infinitely  thin-walled,  perfectly  conducting, 
right  circular  cylinder  Is  shown.  As  Indicated  the  cylinder  is  of 
length  L,  radius  a,  and  height  h above  an  Infinite,  finitely  conducting 
ground  plane.  The  system,  consisting  of  the  cylinder,  ground  plane, 
and  Incident  field  of  electromagnetic  radiation  1s  defined  in  tenns  of 
a combined  cartesian  and  cylindrical  coordinate  system.  The  Incident 
electromagnetic  field  propagates  In  a general  direction  described  by  the 
angle  e with  respect  to  the  z-axis.  The  current  Induced  on  the 
cylinder  by  the  incident  wave  Is  to  be  determined. 

Solution  of  the  above  defined  problem  may  be  greatly  facilitated 
by  the  use  of  Image  and  reflection  theory.  Applying  image  theory,  the 
cylinder,  Incident  TEM  plane  wave,  and  finitely  conducting  ground  plane 
are  to  be  replaced  by  the  object  cylinder,  the  Image  cylinder,  and  two 
free  space  TEM  plane  waves  as  shown  in  Figure  2-2.  Reflection  theory 
can  be  used  to  show  that  the  currents  on  the  Image  are  related  to  the 
currents  on  the  object  by  the  ratio  of  the  complex-reflection  coefficient 
for  the  appropriate  angle  of  Incidence  and  polarization  of  the  Incident 
field  Cl  1-14].  The  reflection  method  cannot  be  expected  to  give  an 
exact  solution  as  the  cylinder  approaches  the  ground  plane  [10], 

[14-16]. 

Two  coordinate  systems  are  defined  by  Figure  2-3,  the  object  and 
the  Image  coordinate  system  Indicated  by  the  subscripts  "O”  and  "1" 
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respectively.  Note  that  since  the  cylinders  are  Identical  and  parallel, 
Zp  and  z^  axes  are  not  necessary.  Let  the  "Incident  field"  be  defined 
as  the  field  plus  Its  reflection  from  the  ground  plane.  The  surface 
currents  Induced  on  the  object  and  Its  Image  are  considered  as  source 
currents  radiating  In  free  space  [7].  The  free  space  Green's  functions 
may  be  used  to  compute  the  scattered  field  at  an  arbitrary  field  point. 

In  terms  of  the  geometry  as  defined  by  Figure  2-3,  appropriate 
magnetic  vector  potentials  for  the  object  and  Image  are: 

Aq(Ro)  ■ magnetic  vector  potential  of  the  object  In  object 
coordinates 

= magnetic  vector  potential  of  the  Image  In  Image 
coordinates 

■ general  field  point  In  space  measured  from  the  object 
coordinate  system 

• same  general  field  point  In  terms  of  the  image 
coordinates 

Ig'  ■ a general  source  point  on  the  object  cylinder  with 
respect  to  Its  coordinate  system 
- a general  source. point  on  the  Image  cylinder  with 
respect  to  Its  coordinate  system 

The  magnetic  vector  potential  of  the  object  Is  expressed  as 

*0 
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Primed  coordinates  Indicate  source  points,  unprimed  coordinates  field 
points; 

space  Green's  function  In  object 
coordinates 

1^{Rq')  " surface  current  density  radiating  In  free  space. 


In  general  the  free  space  Green's  function  has  the  form 


•) 


^-YIFq  - Wq'  I 

l**o  “ ^o'  I 


The  temporal  variation  has  been  suppressed,  where 


(Z-2) 


s - «f  + j«,  (2-3) 

The  complex  frequency  variable,  with 

y • s/c  (2-4) 

c the  speed  of  light  in  free  space. 

By  the  superposition  of  a cylindrical  coordinate  system  on  the 
cartesian  coordinate  system  of  Figure  2-3,  the  following  may  be  shown: 

l>*0  ■ ®o'l  ■ * "o'^  - • ♦o'> 

♦ (2-5) 

where,  Pq'  - a,  a the  radius  of  the  cylinder.  Therefore, 
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,L  2v 


Gq(Pq.4>q.z;Pq'  .♦q*  ,2'  )ad^Q'dz' 


(2-6) 


Reduces  to 


u L 2n  -yR^ 

3^o(Po.*o‘Z)  - ^ 

0 0 ' 


where 


2 . .2 


Rl  • ■ 2p- 


|a  cos(«Q  - *q')  + z - 2')^!^.  (2-8) 


Assuming  no  circumferential  current  components. 


VV'')  ■ 


(2-9) 


Equation  (2-7)  reduces  to 


w I Z-”  ...., 

0 0 ' 


(2-10) 


A similar  process  yields 


p L 2w  -YRg 

(P^.*^.z)  - ^ J f K^(^^.2')^r — ad+;dz'.  (2-11) 

35  ft  "^ft  2 


0 0 


where 
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Rg  ■ - R'JI  ■ [p^  + - 2p^a  cos(«l'^  - ij>j)  + (z  - z')^y 

(2-12) 

The  circumferential  variations  of  the  axial  current  on  an  Infinitely 
long  cylinder  over  a perfectly  conducting  ground  plane  In  a static  mode 
has  been  derived  by  Taylor  [IIQ.  Using  Taylor's  equation. 


. I \2- 


(2-13) 


where 


^0^*6^  " r+'(a/^)cos 


(2-14) 


Iq(2')  ■ axial  variation  of  object  surface  current. 

As  h becomes  large,  the  circumferential  variation  of  the  axial  current 
becomes  uniform  and  thus  (2-14)  reduces  to 


0’*'^  * 2xa 


(2-15) 


Similarly 


1^(2') 


" “Zia 


*2* 


(2-16) 


The  magnetic  vector  potentials  may  now  be  expressed  as 

^ In(2') 

0 0 


w ^ I (2')  “THi 

f f \ir  V 

A n I 
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»^1 »2) 


ad't'^dz' . 


(2-18) 


u ^ I /»•!  “VRfl 

r r *1'’=  ^ e 2 

^ J J RT 

0 0 ^ 


Let  the  general  field  point  of  Figure  2-3  be  located  on  the  object 
cylinder.  With  this  restriction  the  distance  R2  may  be  calculated  In 
terms  of  the  object  coordinate  system;  see  Figure  2-4.  From  the  law 
of  sines, 


sin  ■ a sin  o ■ a sin  ♦q  (2-19) 

and  from  the  law  of  cosines, 

0^  ■ + 4h^  + 4ah  cos  ♦q  (2-20) 

Thus  when  Is  evaluated  on  the  object  surface, 

■ [a^  + 4h^  + 4ah  cos  (2-21) 

and 

♦ ^ - s1n“^  [ ^1  . (2.23) 

L [a^  + 4h^  + 4ah  cos  fg]  J 

Note  when  iK^  Is  evaluated  on  the  surface  of  the  object 

Pq  - a.  (2-23) 

With  the  field  point  on  the  object  cylinder  (2-17)  and  (2-18)  may  be 
written  as 
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ligurtt  2-'4.  Fl«ld  Point  on  th*  Object  Surface 
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I 2ir  , f -yR, 

"1^  / f "lir" 

i\  rk  I 


(2-24) 


0 0 


■ [2a^  - 2a^cos(4Q  - + (r  - z')^] 


(2-25) 


"1^  f f 


*•  ^(z')  ^^'*2 

-J^^~,p-«d4;d2' 

0 0 ^ 


(2-26) 


R2  ■ [2«^  + 4h^  + 4ah  cos  ♦q  - 2[a‘  + 4h‘  + 4ah  cos  ^q]  • 


2 . .u2 


a cos(4^  - ♦^')  + (z  - z')^] 


■ sin 


-1 


a $1n  4, 


[a^  + 4h^  + 4ah  cos  ^q]  J 


(2-27) 

(2-28) 


As  stated  earlier,  the  currents  on  the  object  and  inuge  are  related 
by  the  complex-reflection  coefficient  for  the  appropriate  angle  of 
Incidence  and  polarization  Involved.  If  the  cylinders  are  viewed  as 
consisting  of  a large  number  of  horizontal  dipoles  the  problem  of 
reflection  reduces  to  that  of  a horizontal  dipole  above  a finitely 
conducting  ground  plane.  As  pointed  out  by  Jordon  [14],  In  a plane 
parallel  to  the  axis  of  the  dipole  the  electric  field  Is  parallel  to 
the  plane  of  Incidence,  vertical  polarization.  In  a plane  perpendicular 
to  the  axis  of  the  dipole  the  electric  field  Is  perpendicular  to  the 
plane  of  Incidence,  horizontal  polarization.  The  reflection  coefficients 
for  horizontal  and  vertical  polarization,  as  derived  by  Jordon  D4],  are 
respectively:  13 


(Z-29) 


sin  - [(c  + x)  - cos  ii] 

"h = r-^ 

Sin  i*  + [(e^  + x)  - cos  ti*] 

2 *s 

(e^  + x)sin  ♦ - [(c  + x)  - cos  ♦] 

Ry  - ^ (2-30) 

+ x)sin  ♦ + [(e^  + x)  - cos^  ♦] 


»<here  x-  ~ 

“o 

- relative  permittivity  of  earth 

0 ■ conductivity  of  earth 
Cq  ■ free  space  permittivity. 

Uhen  the  plane  of  incidence  Is  neither  in  a plane  parallel  or  perpendicular 
to  the  axfc  of  the  dipole  a combination  of  vertical  and  horizontal 
reflection  coefficients  must  be  used.  The  reflection  coefficient  for 
this  case  is  given  by 

Rp  ■ cos  0 - R^  sin  0 ...  (2-31) 

The  minus  sign  comes  from  Jordon's  assumed  positive  directions  of 
electric  fields  for  the  Incident  and  reflected  «Mves  [lO.  Figure  2-5 
defines  the  geometry  of  the  reflection  coefficient.  Since  uniform 
circumferential  current  has  been  assumed  no  generality  Mill  be  lost 
hy  restricting  the  field  point  of  the  object  cylinder  to  lie  along 
the  line  fg  ••  0.  With  this  objective  and  image  current  related  by 
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Source 


Field 


of  Reflection  Coefficient 


S 


i^{z')  - Rpyz') 


(2-32) 


It  win  be  necessary  to  express  the  reflection  coefficient  In  terms  of 
the  object  coordinate  system.  With  this  objective  In  mind  the  angle 
and  0 defined  In  Figure  2-5  can  be  expressed  as  functions  of  and 
z such  that, 


♦ ■ tan 


-1 


2h  + a(l  + cos  ♦q') 

L sln^  + (z  - z')^)^ J 


(2-33) 


and 


0 ■ tan 


-1 


4 • 
0 


a sin  ^ 


(2-34) 


Substituting  for  and  0 In  equations  (2-29),  (2-30),  and  (2-31) 
yields 


1 - 


1 + 


(ty. 

sln^  » 


+ X - 1)  1 ** 

+ 1 J 

sin^  ^ 


a sin 

(a^  sln^  + (z  - z')^) 


(2-36) 


{‘y.  + X) 


L sift^  » J 


jjLiUJl 


(c  ,x)  i1  (aWfiMz-z')2) 

sln^  ♦ J 


Dropping  the  current  subscripts  of  equation  (2-32)  the  magnetic  vector 
potentials  become 


adiji^clz' 


(2-36) 


Ao2(4o*^) 


»^0 


2it 


f 


RpKz')  c‘^^2 
2«a 


ad^jdz' 


(2-37) 


The  total  field  consists  of  the  Incident  and  scattered  field, 
where  the  scatteied  field  Is  that  due  to  the  current,  I(z').  To  Insure 
uniqueness  the  total  field  must  satisfy  certain  boundary  conditions 
D&]. 


‘ fine  * f,- 


(2-38) 


The  scattered  field  Is  related  to  the  Induced  currents  and  charges  by 


ts(R)  - -sA^jdO  - Wj(ir),  (2-39) 

where  “ the  total  "scattered"  magnetic  vector  potential 

♦j.  ■ the  total  "scattered"  electric  scalar  potential, 


which  can  be  related  to  through  the  Lorentz  gauge  condition, 

. V . HAJO 

♦‘(f)  ■ -Tcnri-  • (2-'"» 


*0‘0' 


Equation  (2-41)  becomes 


t_(lf)  - -s  r v[v  . HAm] 

L W J 


(2-41) 
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since  y • ^ and  c - I/^WqCq,  y » s UqCq.  Then, 

r,(1f)  - -s  ^ v[7  . . (2-42) 

On  the  surface  of  the  object  cylinder,  the  boundary  condition  Is 

n X F-f  - 0,  (2-43) 

with  n being  the  outward  normal  unit  vector  on  the  cylinder  surface. 

An  equivalent  representation  of  the  boundary  condition  Is 


’Inc 


tan 


1 "S.nl 


(2-44) 


Which  means  that  the  tangential  components  of  the  incident  and  scattered 
fields  must  cancel  on  the  object  surface  for  boundary  conditions  to  be 
satisfied.  With  the  above  consideration  equation  (2-42)  becomes 

"'"'UnJ  . (2-45) 

*0 

since 


5^j,(Fq)  ■ ^1  " ^0  ^1  *z* 

L 2x  -yR 

‘ If  nq” 


(2-46) 


L 2* 


-yR^ 


adfjdz' 


(2-47) 


0 0 
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Substituting  (2-47)  into  (2-45)  one  arrives  at 


"inc 


tan 


I • s[^l  j 5? 


(2-«) 


L 2if  -YR 


1 


\j  f ^ her  * j j -wi — 

loo  ‘ 00  ^ 


> p RA(z') 

J J -Ai— ^r-acl*)dz' 


Since 


[l  - ^v(7-)]  • -»*[l  - ^7(7.)] 

[7(7.)  -T^]  (2-49) 


-S^ 

■? 


[2-48]  becomes 


-4weoS 


J ) • 


L 2if 


[v(V7)  - Y^] 


[I  I % 

' 0 0 


-yR, 


=-1  . 

ira  Ri 


ad*^(li'  + f / 

•'O  *'0 


Rpl(2')  /^'^2 

•~2ii 


idi'J  , 


(2-50) 


The  differential  operator,  [^(7*}  - Y^]i  In  the  above  equation 
reduces  to 


(&■■■) 


(2-61) 
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since  equation  (2-50)  has  only  a 2 component.  Finally  (2-50)  reduces 
to 


k HMz-)  -’"2 

* J J “25i RT"'** 

*'0  ^ 


;d2'|  . 


(2“52) 


The  Incident  field  as  shown  In  Figure  2-2  can  be  expressed  as 

A A 

Eq  exp[-Y(z  cos  0 - X sin  0)][cos  ® ^ ® 

- exp[-Y(z  cos  0 + X sin  o)](cos  ^ - sin  0 a^,) 

(2-53) 

where  Is  as  defined  In  equation  (2-30).  The  minus  sign  again  comes 
from  assumed  positive  directions  of  electric  fields  for  the  Incident 
and  reflected  waves  In  the  derivation  of  R^  D4].  Since  only  the  axial 
component  of  the  Incident  field  Is  to  be  used, 

E ■ Eq  exp[-Y(2  cos  6 - X sin  e]s1n  0 
*1nc  ^ 

+ RyE^  exp[-Y(z  cos  0 + x sin  0)]s1n  0.  (2-54) 


The  axial  component  of  the  transmitted  field  Is 
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(2-55) 


E,  ■ (1  - Ru)Ert  exp[-Y(z  cos  0 - X sin  0)3s1n  0 
^tran  ^ ^ 


To  determine  the  fields  must  satisfy  the  boundary  condition 


^tan 

*t/an 


X 


-h 


(2-56) 


The  Incident  field  was  defined  to  contain  the  direct  and  reflected 
field,  therefore 


(1  - Ry)EQ  exp[-Y(z  cos  0 + h sin  0)]s1n  e - 


Eq  exp[-Y(z  cos  0 + h sin  6)]s1n  0 + 


RyE^  sin  0 exp[-Y(z  cos  0 - h sin  o)]  (2-57) 


And  upon  simplification 


E^  ■ -Eq  exp[-Y  2h  sin  0].  (2-58) 


Substituting  (2-58)  into  (2-54),  and  evaluating  E_  on  the  cylinder 

^Inc 

axis  (approximate  kernel),  Instead  of  on  Its  surface. 


E^  ■ Eq  exp[-Y  z cos  0]s1n  0^1  - R^  exp[-Y  2h  sin  0]  j . 

(2-B9) 


As  stated  earlier  we  have  restricted  the  field  point  on  the 
object  to  lie  along  the  line  ■ 0.  With  this  restriction 

R-j  ■ [2a^  - 2a^  cosC+q  - fQ*)  + (z  - z')^]  (2-60) 
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and 


where 


Reduces 


Let 


and 


$0 


*4 

Rg  “ [2a^  + 4h^  + 4ah  cos  - 2(a^  + 4h^  + 4ah  cos  ^Iq) 


a cos(4^  “ ♦!)  + (z  - z‘ )^] 


(2-61 ) 


♦ ^ « sin 


•1  ® ♦o 

. [a^  + 4h^  + 4ah 


cos  J 


(2-62) 


to 

'*1  “ '^1 

^Vo-0“ 


rg  - R^ 


- [28^(1  - cos  I-')  + (z  - z')^]' 


■ [(2a^  + 4ah){l  - cos  ^\)  + 4h^  + 

•A 


(2-63) 


■♦o-o“ 


9 ^ 

(2  - Z')^]  . 


(2-64) 


2 

1 - cos  « 2 sln^  -f 

(2-65) 

d diameter  of  the  cylinder 

- 2a 

(2  66) 

r,  . [i,2  sln^  (2  - i-)2] 

(2-67) 

4 ' ^ 

rg  ■ [(d^  + 8ah)s1n^  -^  + 4h^  + (z  - z' )^] 

(2-68) 

Also  defining  two  functions 
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"vr. 


r e ^ 

fo(z,2',s)  - J -f- 

0 r, 


(2-69) 


and 

f^(z,2'.s)  - 
Equation  (2-52)  becomes 


/Cl 


'r,  ad4*i 

0 rg  1 


(2-70) 


(2-71) 


Where  E Is  given  by  equation  (2-59).  In  order  to  better  represent 
*1nc 

the  complex  frequency  dependence. 


L 

f Cfo(z.2'.s)  + RFfi(z.z’.s)]dz‘. 

(2-72) 

This  Is  the  Integro-differentlal  equation  to  be  solved  for  the  unknown 
Induced  current. 
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Application  of  the  Moment  Method 
The  Integro-differentlal  equation  shall  be  cast  Into  matrix  form 
suitable  for  a numerical  solution;  this  process  Is  called  the  method  of 
mements  DS,  6cQ,  Bl  -22].  The  portion  of  the  cylinder  at  ♦q»0,  “a  thin 
wire".  Is  broken  Into  N small  zones,  with  the  Induced  current  on  each 
of  these  zones  assumed  to  be  constant.  The  Integration  over  the  "wire" 
Is  now  (ipproxlmated  by  the  sum  of  Integrals  over  N segments. 

Current  on  the  wire  may  be  described  by  the  function 


where 


and 


l(z',s)  - I %(s)In(z') 
n 


(2-73) 


a m unknown  coefficient  of  constant  current  on  n^^  subsection 
n 


1 , for  z”  < z < z”"^^ 
0,  elsewhere 


(2-74) 


In  Figure  2-6  Is  shown  the  geometry  pertaining  to  the  method  of 
moments.  The  boundary  condition  that  current  at  wire  ends  be  zero  Is 
automatically  satisfied  by  allowing  the  two  end  zones  to  extend  past 
the  surface  of  the  cylinder.  Zone  length,  match  points,  and  end  points 
are  given  respectively  by 

A ■ U/(N-1)  ■ length  of  a zone  (2-75) 

z,.  " (m-l)A,  m ■ 1,2,...,N 
z"  - (n-3/2)a,  n - 1,2,...,N+1 


24 


(2-76) 

(2-77) 


BirfiSECTION  Eims 

X M/VTCH  POIIITS 

1,1  « (m~l)A  (HATCH  POIIITS) 

B " .N 

(11  . (n..3/l*)A  (8UBSECT101I  ENDS) 
n 

a " L/(M«-1)  LENGTH  OF  ZONE 
H • HUMBER  OF  SUBSECTIONS 


J 

■ « 

• 

J 

;»2 

*1-0 

Fieurc  2-6  . ttomont  Method  Partitioning  of  Qoometry 
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Mhere 


L ■ length  of  wire 
N ■ number  of  zones. 


With  current  on  the  wire  expressed  as  a pulse  function,  the  integro- 
dlfferentlal  equation  (2-72)  may  be  written  as 

(-AxcqS)  jEj,  sin  e e"^«°*®Cl-RXY2hs1ne3  j . 

^n+1 

* (2-78) 

This  equation  Is  to  be  satisfied  at  discrete  match  points,  these  points 
being  described  by  equation  (2-76).  Derivatives  In  equation  (2-78)  win 
be  approximated  by  a finite  difference  technique,  that  is 


[F(z  + m)  - 2F(z)  + F(z  - Az)]. 

dz^  (Az)"^ 


(2-79) 


Using  (2-79)  In  (2-78)  yields 

V>  !'o 


(-4...S)  Ie.  stn  e [i  . r„ 


.n+1 


n A %n  « 


♦ Fo(*,„.l »*'.«)  - RfrCF^(z^,,r'.s)-(vV+2)F^(*^.*',s) 


n - 2,3 N-1 


♦ F4(z„.it*'i»)]  dz' 

’ ” ' M ■ 2,3, . . . ,N“1 


26 


(2-80) 


where 


1 . l] 

a sin  1>q' 

1 ♦ [14jU.  i]'’ 

A*[i4jli+  l]"* 

[a2s1n^^Q'+(2^-z'] 

* ■ V * -k 


B - -tan"’ 


2h+a(l+cos^Q' ) 
[a2stn%'+(z^-z')2] 


This  equation  may  be  manipulated  Into  the  form 
\r(s)  - f(s)T(s) 


(2-81) 


(2-82) 

(2-83) 


(2-84) 


where  a single  bar  represents  a column  matrix  or  vector  and  a double 
bar  Indicates  a square  matrix. 

Let  V(s)  ■ the  source  vector  ■ 

where  v„  - the  matrix  elements  of  ^(s) 

m 

- (-4*Cqs)  j Eq  sin  e 

m-2.3 N-li  (2-85) 

7(s)  - the  Impedance  matrix  - 
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where 


*»».  “ the  matrix  elements  of  Z(s) 
wn 

-n+1 


jfn  ^ 


F,(V'''S)+fl(Vl>^’-*>3|'''' 

n ■ 2,3,. . . ,N-1 
m • 2,3,. .. ,N-1 


(2-86) 


and  T(s)  ■ the  response  vector  - [1^], 

where  ij,  ■ the  matrix  elements  of  T(s) 

■ Ojj,  unknown  coefficient  of  constant  current  in  the 

n^^  zone  n ■ 2,3,..., N-1  (2-87) 


The  unwieldy  appearance  of  equation  (2-86)  may  be  somewhat  improved  by 
defining  two  functions 

.n+l 


with  r^  ■ [d^  $in^  ♦ (z^  - z')^3 


and  H 


z"*^  2*  . 

r c 


(2-88) 


(2-89) 


(2-90) 
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with 


(2-91) 


4 ' lg| 

Tg  ■ [(d^  + 8ah)s1n2  -J-  + 4h^  + (z^  - z')^]  . 

Now  equation  (2-86)  may  be  redefined  as 

c"a  ( ffl  n m n 

ii  n Ri  n m n 

+ H.  (2-92) 

m n ) 

When  ♦q'^O®  and  z'“z^  then  r.j-0,  (see  equation  (2-89))  thus  the 

Integrand  of  Hq  equation  (2-88),  will  be  singular.  Tesche 

m n 

C!3  ] provides  methods  of  evaluating  Integrals  of  this  type.  Note  that 

2 

the  Integrand  of  (z„is)  Is  never  singular  due  to  the  term  4h  In 
^m  n ™ 

equation  (2-91). 

At  the  Integrand  of  Hq  (Zuj^.|,s)  is  singular  when  Zj^^^-z', 

n n+1  ^ ^ 

where  z ^z'  , Hence,  at  this  point. 


* - *mfl  - ^ 


(2-93) 


Using  equation  (2-76)  and  (2-77)  It  Is  seen  that  the  singularity  occurs 
when 


-3/2  < (m-n)  < -1/2.  (2-94) 


M and  n are  integers,  therefore  Hq  (z^is)  Is  singular  when  m-n  * -1  or 

ro  n 

when  m>n~1.  A similar  argument  may  be  applied  to  show  that  Hq  (Zm>s) 

''ran"' 

is  singular  when  m-n  and  that  Hq  (z  ^^,s)  Is  singular  when  mon-M. 

m n 
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Following  the  procedure  outlined  by  Tesche  [23],  let  the  Integral  of 

Hq  singularity  be  T, . 

*1 

Ti  ■ Hq  (2-96) 

n 

where  Kq  (2^*^)  9lven  by  equation  (2-88).  Transform  variables 
m n 

In  (2-88)  as  follows: 
let  ♦ - ♦q' 

z ■ - z'  (2-96) 

dz  ■ -dz' 


and  the  limits  of  Integration  become 


m 4 

in*n  m*n 

- &/2 

1 . - z"^l 

Jm-n  •" 

1 - -a/2 

'uwn 

(2-97) 

Thus, 

- 7 0 

di^dz 

(2-98) 

where 

r^  « [z^  ♦ d^  sln^  e/2]** 

Define 

h ■ 

when  m > n-1 

(2-99) 

and 

T3-\„<Vr*) 

when  n • n+1 

(2-100) 

it  can  be  shown  that  T|  "12"  T^. 
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Equation  (2-98)  may  be  prepared  for  numerical  evaluation  by  following 

a procedure  similar  to  that  used  by  Tesche  [23].  A Maclaurln  series 
-yr, 

for  e Is  as  follows: 


-yr. 


2 2 3 3 

r,  Y r,  Y 


. 1 . , ....  '1  ' . 
' -l-vr,  ... 


Using  only  the  first  two  terms  of  (2-101)  equation  (2-98)  may  be 
written  as 


(2-101) 


Tl  - 


(2-102) 


'0  Si/ 2 


also 


A/2  ^ 


(2-103) 


and  therefore 


■ 2 / *"j  JI*  [ I]  j * 


2ii 

2 in(s1n  |)d*^  - 


2»Ay. 


(2-104) 
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The  first  term  In  the  above  equation  Is  non-singular  and  can  be  evaluated 
numerically,  the  second  term  Is  easily  evaluated  analytically  be  expressing 
%n(s1n  t/2)  in  a series,  that  Is. 


in{s1n  ■ -in  2 - £ ^ cos  n* 

n-1  " 


so 


Finally. 


i 


in(s1n  4/2) 


- 1 /" 

■ -2*  in  2 - E ^ / 

n-1  " 1 


(2-105) 


cos  n4d4  (2-106) 


T,  - 2 


I / 


+ sin  4/2  1 d4 


I'} 


+ Tt(in  4 - ay)  . 


(2-107) 


The  Iritegro-differentlal  equation  has  been  cast  as  a system  of 
matrices,  which  may  be  solved  for  the  unknown  Induced  current. 


Application  of  the  Singularity  Expansion  Method 
The  singularity  expansion  method  (SEM).  formalized  by  Baum  [23-25], 
and  applied  by  many  others  [1-3],  [26].  [27-30],  Is  a method  for  charac- 
terizing the  response  of  scattering  objects  Illuminated  by  either 
transient  or  steady*state  electromagnetic  radiation.  A general  outline 
describing  the  theory  of  SEM  will  be  presented  with  special  attention 
given  to  areas  of  Interest. 
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The  time  domain  current  may  be  determined  from  a knowledge  of  the 
exterior  natural  resonant  frequencies  of  the  scatterer  In  much  the  same 
manner  as  In  classical  circuit  theory.  The  natural  frequencies  of  the 
“wire,"  Sg,  are  those  such  that  the  homogeneous  version  of  equation 
(2-84) 

Z(s)T{s)  - 0 (2-108) 

has  a non-trivlal  solution  for  I(s^jj).  This  Implies  that  the  determinant 
of  ? must  vanish  at  these  complex  natural  frequencies;  and  therefore  the 
equation  for  determining  these  resonances  becomes 

det  Z(Sj,)  - 0 (2-109) 

The  natural  frequencies  must  lie  1n  the  left-hand  portion  of  the 
s-plane  since  exponentially  Increasing  currents  are  not  allowed.  Further, 
the  poles  must  occur  in  conjugate  pairs  since  the  time  domain  current  is 

real.  No  poles  may  reside  on  the  jw  axis  since  the  radiation  process 

requires  that  the  current  on  the  wire  be  eventually  rer-o.  It  1s  assumed, 
without  proof,  that  the  poles  are  all  simple.  This  has  been  substantiated 
numerically  [30]. 

The  matrix  equation, 

. 7(5)T(s)  - V(s)  (2-110) 

has  as  a solution 

T(8)  ■ y ■‘'(s)\r(s).  (2-111) 
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The  time  domain  response  Is  given  by  the  Inverse  laplace  transform 
In  the  form 


Hi) 


1 


Z "^(s)7(s)e*^ds. 


(2-112) 


Instead  of  numerically  evaluating  this  Integral  along  the  Oq 
contour,  SEM  assumes  that  the  Inverse  of  the  Impedance  matrix  may  be 
expressed  as 


? '^(s)  ■ I . (2-113) 

This  assumption  allows  evaluation  of  equation  (2-112)  by  the  Cauchy 
Residue  Theorem.  Equation  (2-113)  Is  a summation  over  all  poles  In  the 
complex  s-plane.  R^j  Is  defined  as  the  system  residue  matrix  at  the  polo 
s^jj.  This  residue  matrix  Is  a dyadic  and  can  be  represented  as  the  outer 
product  of  vectors  Independent  of  s as 

R^  • (2-114) 

where  Is  the  natural  mode  vector  and  Is  a solution  to 

Z(s„)R^  • 0 (2-115) 

and  Cy  Is  the  coupling  vector,  which  satisfies  the  equation 

- 0 (2-n6) 
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(T  denotes  transpose).  For  the  electric  field  formulation,  where  the 

mm 

Impedance  matrix  Is  symmetric,  i.e.  (Z(s^)  " Z (s^j^)),  the  natural  mode 
vector  and  coupling  vector  are  equal,  and  therefore 

(2-117) 

Let  (T^  be  normalized  such  that  its  maximum  element  1s  real  and  equal  to 
unity. 


[Ca  ] 

normalized 


(2-118) 


and  define 


normalized 


(2-119) 


such  that 


(2-120) 


where  Bj,  Is  the  normalization  coefficient.  Since 


Z "^s) 


Bir  r ' 

° “o 
s-s„ 


(2-121) 


The  normalization  coefficient  Is  determined  through  consideration  of  a 
particular  singularity  [2],  Sp,  such  that 

Cp''‘?(s )?  “ ’ (s  )fp  - fp^U  Cp  - Ifp^e'p  (2-1 22 ) 

(NOTE:  U Is  the  Identity  matrix.) 
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Substituting  equation  (2-121)  Into  (2-122)  yields 


f /2(s)  I 

K o a 


C 

P P 


or 


F '^Z(s)f  C 

p * * o g p 

s-s 
* o 


C . 
P P 


(2-123) 


Note  that 


so  (2-123)  may  now  be  written  as 

; p ■ g V 

J “o  S-S^  P P 


(2-124) 


(2-125) 


By  definition, 


Z(s)  - ?(Sp) 


(2-126) 


Where  V (Sp)  ■ 

S-Sp 

And  6 . is  the  Kronecker  delta  function.  Taking  the  limit  as  s approaches 
pa 

Sp  In  (2-125)  results  In 


Bf^Z'Cs  )C’C’^7 
Vp  ' P^  P P P 


CjlT. 


P* 
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(2-127) 


From  whltii  we  conclude  that 


Finally, 


p r (s  Yr  * 

p ' p'  p 


VJ" 


1 f 

■ 5iT  / I ^ 

On-J- 


where  N„  ■ 


(2-128) 


(2-129) 


Is  defined  as  the  coupling  coefficient.  With  the  matrix  elements  of 
V(s)  defined  by  (2-85),  let  the  Incident  field  be  a step-function  plane 
wave,  such  that 


Eo(s)  • Eg/S 

The  matrix  elements  of  V(s)  are  then  defined  as 


Vm  - -4’'Cg  jEg  s1n  e e [1  - ®]  j (2-130) 

where  m > 2,3,. . . , N-1 , 

Evaluating  (2-129)  through  the  residue  theorem  will  produce  appropriate 
Heaviside  functions,  which  are  viewed  as  enforcing  causality  [30],  [1-3], 
The  exponential  dependence  of  (2-129)  Is  expressed  by 

(2-131) 
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where 


D ■ -A»‘iQEQ  sin  0. 


Equation  (2-131)  can  be  written  as 


V ■ D 
m 


s[t  - z /c]  s 
s "*  -R^e 


And  now  (2-129)  becomes 


[t  -^2 


+2h  sin  e 


i(t) 


2iJ 


/. 


® “n  “n 


7^(5)  Is  a vector  with  matrix  elements  defined  by  (2-132). 
simply, 

V„(s)  - [v„e*^]  - 7(s)e*^ 

If  we  let 


'1  ■ '111'= 


and 


'2- 


z_  + 2h  sin  0 
m 


Then  the  matrix  elements  of  \{i)  are 


v_e*^  ■ D[e 
m 


s(t--«^)  S(t-T,) 


Ry® 


3. 


(2-132) 

(2-133) 

Stated  more 

(2-134) 

(2-135) 

(2-136) 

(2-137) 


Evaluation  of  (2-133)  may  be  carried  out  by  using  the  residue  theorem. 
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Note  that  the  complex  natural  frequencies*  natural  mode  vectors, 
and  normalization  coefficients  are  not  functions  of  the  Incident  field. 
Only  \{t)  Is  altered  upon  a change  In  the  angles  of  Incidence.  Conse- 
quently, once  s , 0 , and  C.  are  found  for  a particular  L,  h,  and  a, 

u a 

the  current  excited  by  any  Incident  field  Is  easily  found. 

The  singularity  expansion  method,  applied  to  this  transient  electro- 
magnetic problem,  produces  a system  of  matrices,  which  are  solved  for 
the  Induced  current  on  the  object  cylinder  as  a function  of  time. 

Approximations  and  Limitations 

A summary  of  the  approximations  used  thus  far  will  be  presented 
so  that  the  results  obtained  may  be  viewed  In  the  proper  perspective. 

The  assumptions  and  approximations  are: 

1.  Current  Is  assumed  to  flow  only  In  the  direction  of  the  cylinder 
axis. 

2.  Boundary  conditions  are  enforced  only  on  the  axial  component 
of  the  tangential  electric  field. 

3.  End  caps  on  the  cylinder  are  Ignored. 

4.  Current  Is  assumed  to  be  uniformly  circumferentially  distributed. 

5.  The  moment  method  Is  an  approximate  numerical  solution. 

6.  The  "T"  function,  equation  (2-107),  for  evaluating  the  singular 
Integral  Is  an  approximate  solution. 

7.  The  reflection  method  does  not  account  for  the  currents  Induced 
in  the  earth. 

Assumptions  one  and  two  require  that  the  cylinder  be  thin,  L » a. 
The  third  assumption  Is  valid  If  a « this  Is  seen  to  be,  because 
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current  Induced  on  the  end  caps  will  not  significantly  contribute  to  the 
scattered  field  [31].  Uniform  circumferential  variation  of  the  current 
Is  a valid  approximation  when  the  cylinder  Is  many  radii  away  from  the 
ground  plane.  The  assumptions  involved  In  Tesche's  solution  for  the  T 
function  places  two  restrictions  on  the  problem  [2|.  The  main  restriction 
requiring  L » a,  and  secondly  this  formulation  Is  not  applicable  to 
high  frequency  analysis. 

The  reflection  method,  although  valid  when  the  cylinder  Is  located 
far  from  the  ground  plane,  Is  not  valid  when  the  cylinder  Is  close  to 
the  ground  plane.  Additional  terms  must  be  added  to  the  equations  to 
account  for  the  "surface  wave"  [14].  The  term  surface  wave  as  used  here 
Is  as  defined  by  Norton  [32].  In  a dissertation  by  Jerry  McCannon  p6], 
a study  was  done  of  a vertical  dipole  using  several  different  formulations. 
The  reflection  method  was  found  to  give  answers  within  1 to  2 percent 
of  the  exact  solution  when  the  height  of  the  dipole  was  greater  than 
3/8  X. 

In  a recent  Investigation  by  Sarkar  and  Strait  D0]i  it  was  found 
that  the  reflection  method  for  the  horizontal  electric  dipole  produced 
results  within  lOX  of  the  exact  Sommerfeld  formulation  If 

h > X,  (2-138) 

>Fr 

where 

h •>  height  of  the  dipole  above  the  ground  plane 
■ relative  permittivity  of  ground  plane 

X > free  space  wave  length. 
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Throughout  this  study  an  average  of  15  was  used  for  relative 
permittivity  of  earth  In  this  case,  which  results  In  h > .065  X. 


III.  NUMERICAL  RESULTS 


A computer  code  has  been  written  to  implement  the  equations  developed 
In  the  previous  section.  This  code  Is  used  to  determine  the  natural 
resonances  and  natural  mode  vectors  of  the  scatterer.  The  scatterer  Is 
described  by  a general  length,  radius,  and  height  above  the  ground 
plane.  The  finitely  conducting  ground  plane,  assumed  nonmagnetic  v * Wq, 

Is  characterized  by  o,  the  conductivity,  and  c ■ the  permittivity. 
Singularities,  1.e.( natural  resonances),  occur  In  layers  In  the  complex 
plane  an1  will  be  described  by  s.  „ where  denotes  the  layer  of  the 
pole  and  “n"  the  pole  within  the  layer,  this  method  of  description  Is 
that  used  by  Tesche  [30].  This  data  will  present  trajectories  of  the 
first  three  poles  of  the  first  layer,  s^^,  s^g,  and  s^j. 

Figure  3.1  shows  the  movement  of  the  singularity  s^.|  In  the  complex 
plane  as  several  system  parameters  are  varied.  The  outer  dashed  spiral 
through  the  points  labeled  A describes  the  movement  of  the  pole  s^-j  as 
the  scatterer  Is  brought  near  the  ground  plane,  the  conductivity  Is  held 
constant  at  o > 120.0  along  this  curve.  It  may  be  concluded  that  the  Q 
of  the  scatterer  stabilizes  as  h/i  becomes  large  since  the  dashed  spiral 
seams  to  become  confined  to  a smaller  area  In  the  s-plene.  For  each 
value  of  h/t,  the  conductivity  Is  varied  from  cr  • 120,  point  A on  the 
dashed  curve,  to  o ■ 1.2  x 10'^,  point  Qi  Intermediate  values  are  shown 
at  points  B through  F.  The  value  o • 1.2  x 10"^  and  o « 1.2  x 10*^ 

at  points  E and  F respectively  correspond  to  typical  values  of  condKtIvlty 
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for  normal  terrain.  Let  the  paths  traversed  by  the  pole  for  a given 
value  h/t  be  called  the  Inner  spirals.  It  Is  seen  that  each  Inner 
spiral,  corresponding  to  a given  value  of  h/t,  converges  to  point  G 
as  the  conductivity  of  the  ground  plane  1$  reduced,  this  result  Is  to  be 
expected  since  point  G Is  the  location  of  s^^  for  free  space  conditions 
[29].  When  the  pole  Is  displaced  from  position  G along  one  of  the  Inner 
spirals  more  energy  Is  being  reflected  from  the  ground  plane,  this 
follows  front  the  fact  the  Incident  Illumination  experiences  a greater 
discontinuity  as  the  value  of  a Is  Increased.  In  the  limit  as  a becomes 
very  large  all  the  Incident  energy  would  be  reflected  and  the  problem 
becomes  that  of  a cylindrical  scatter  over  a perfectly  conducting 
ground  plane.  Figure  3-2  Is  also  a plot  of  as  the  value  of  h/t 
and  a vary,  the  relative  permittivity  Is  held  at  five.  Unlike  Figure 
3-1,  point  G Is  not  the  same  In  the  limiting  case  of  small  o,  but 
rather  each  Inner  spiral  converges  to  some  point  along  an  inner  dashed 
spiral.  Although  the  conductivity  becomes  small  the  relative  permittivity 
remains  at  five,  thus  the  Incident  wave  will  meet  some  discontinuity 
and  there  will  be  energy  reflected  from  the  ground  plane.  Figures  3-3 
and  3*4  display  similar  Information  as  described  above.  In  Figure  3 -3 
and  3.4  the  relative  permittivity  Is  fifteen  and  one  hundred  respectively, 
as  0 again  gets  small  point  B Is  further  from  that  of  the  free  space 
position  due  to  the  larger  vilues  of  To  sunnarlze,  displacement 
of  the  singularity  from  its  free  space  position  Is  a function  of  the 
magnitude  of  discontinuity  In  the  ground  plane  whether  It  be  produced 
by  the  conductivity  or  relative  permittivity.  Figures  3-5  through  3-8 
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may  be  interpreted  similarly.  Note  that  In  Figures  3-5  and  3-7  the 
Inner  spirals  cross  the  outer  dashed  spiral,  this  Is  most  probably  due 
to  numerical  error.  When  the  method  of  moments  Is  employed,  a rule  of 
thumb  Is  to  choose  the  number  of  lones,  N.  In  accordance  with 
N “ 10  X sL/^rc.  This  Implies  that  if  s ■ ju,  ten  cells  per  half- 
wavelength are  used  [30].  This  rule  was  adhered  to  for  pole  s,^,  the 
first  resonance  of  the  scatterer,  but  N was  not  changed  for  s-j^,  or 
s^^  this  means  that  there  were  only  5 zones  per  half  wavelength  for 
s^2*  '^nd  3 zones  per  half  wavelength  for  s^^. 

The  real  and  Imaginary  part  of  the  normalized  mode  vector  for  the 
singularities  S|i,  and  s^^*  shown  In  Figures  3-9  through  3-12.  Note 
that  the  modes  are  either  even  or  odd  functions  about  the  scatterer 
midpoint.  It  is  also  seen  that  the  Imaginary  part  of  the  mode  vector 
Is  much  less  than  the  real  part.  Indicating  the  modes  are  almost  real 
functions  of  position.  The  conductivity  and  relative  permittivity  were 
also  varied  when  calculating  the  normal  mode  vectors,  results  showed  the 
normal  mode  vectors  being  relatively  Independent  of  these  quantities. 
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